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Abstract 

Aluminum hydride (alane) AIH3 is an important material in hydrogen storage applications. It 
is known that AIH3 exists in multiply forms of polymorphs, where a-AlHa is found to be the most 
stable with a hexagonal structure. Recent experimental studies on 7-AIH3 reported an orthorhom- 
bic structure with a unique double-bridge bond between certain Al and H atoms. This was not 
found in a-AlH3 or other polymorphs. Using density functional theory, we have investigated the 
energetics, and the structural, electronic, and phonon vibrational properties for the newly reported 
7-AIH3 structure. The current calculation concludes that 7-AIH3 is less stable than a-AlH3 by 2.1 
KJ/mol. Interesting binding features associated with the unique geometry of 7-AIH3 are discussed 
from the calculated electronic properties and phonon vibrational modes. The binding of H-s with 
higher energy Al-p, d orbitals is enhanced within the double-bridge arrangement, giving rise to a 
higher electronic energy for the system. Distinguishable new features in the vibrational spectrum 
of 7-AIH3 were attributed to the double-bridge and hexagonal-ring structures. 

PACS numbers: 71.20.Ps, 63.20.Dj, 61.66.Fn, 61.50.Lt 
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I. INTRODUCTION 



Aluminum hydride AIH3 (alane) is among the most promising metal hydrides for the hy- 
drogen storage medium, containing a usable hydrogen fraction of 10.1 wt.% with a density 
of 1.48 g/ml. The enthalpy of the reaction is found to be low, resulting in minimal heat 
exchange for both charging and discharging reactions. This material is thermodynamically 
unstable near ambient conditions, but it is kinetically stable without releasing much hydro- 
gen for years. However, extremely high hydrogen pressure (exceeding 25 kbar) is required 
to achieve charging. The decomposition of AIH3 occurs in a single step: 

AlH, ^Al + ^H2. (1) 

How to overcome the kinetic barriers and to find new routes to synthesize AIH3 for reversibil- 
ity have been the focus of many research activities recently. 

Early studies identified seven polymorphs of AIH3: a, a', (3, S, e, 7, and It was 
suggested experimentally that the a phase is the most stable, followed by the (3 and 7 
phases.-i^ The decomposition of the 7 and f3 phases is faster than that of a phase. An 
reaction enthalpy of 7.1 KJ/mol-H2 and 11.4 KJ/mol-H2 for 7 and a, respectively, is 
reported.-*^ The exothermic transition to the a phase from the 7 phase is also reported. 
Only the a polymorph AIH3 has been experimentally investigated extensively, including 
structural characterization,- thermodynamic measurements,-!^ and thermal and photolytic 
kinetic studies .-'^'^'^^lii'i^ Theoretical calculations^^ii^ii^ of a-AlH3 from first principles have 
been performed to study the structural stability and electronic and thermodynamic proper- 
ties. Limited studies have been conducted on other polymorphs and their properties. 

Recently, the crystal structure of 7-AIH3 is reported by two separate groups using syn- 
chrotron X-ray powder diffraction^ and powder neutron diffraction,— respectively. A unique 
feature of double-bridge bonds involving A1-2H-A1 is identified in addition to the normal 
bridge bond of Al-H-Al as found in a-AlH3. In the present study, using density functional 
theory and the linear response approach, we investigate the electronic properties and phonon 
spectra for the newly published 7-AIH3 crystal structure.—*^ The phase stability and inter- 
esting binding characteristics in 7-AIH3 are presented and compared with that of a-AlH3. 
The origin of the cohesive energy difference in these two phases is discussed. Distinct phonon 
vibrational modes arising from the double-bridge and hexagonal ring structures are identi- 
fied. 
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II. CALCULATIONAL PROCEDURES 



The calculations are based on density functional theory.— The Kohn-Sham equations are 
solved in a plane- wave basis using the Vienna ab initial simulation package (VASP).— 
For the exchange-correlation functional, the generalized gradient approximation (GGA) of 
Perdew and Wang (PW91)^ is employed. The electron-ion interaction is described by 
ultrasoft pseudopotentials (USPP).— The k space integrals are evaluated using the sampling 
generated by the Monkhorst-Pack procedure.— The calculation for both a- and 7-AIH3 
structures are performed with a A;-point mesh of 7 x 7 x 7. The relaxations of cell geometry 
and atomic positions are carried out using a conjugate gradient algorithm until the Hellman- 
Feynman force on each of the unconstrained atoms is less than 0.01 eV/A. The nuclear 
coordinates are first allowed to relax while the cell volume is fixed at the experimental 
value. Then simultaneous relaxations of the cell volume, shape and atom coordinates are 
conducted. For the total energy calculations, the plane-wave energy cutoff is 600 eV. The 
self-consistent total energy converges to within 10~^ eV/cell. 

The vibrational properties are studied with density functional perturbation theory within 
the linear response.— The dynamical matrices are obtained for a uniform grid of q vectors 
of 4 X 4 X 4 and 3x3x3 over the Brillouin zone of a- and 7-AIH3, respectively. This 
dynamical matrix is then Fourier-transformed to real space and the force- const ant matrices 
are constructed, which are used to obtain the phonon frequencies. The PWSCF numerical 
code^ was used in our calculations for the zero-point energies and phonon density of states. 

III. RESULTS AND DISCUSSIONS 
A. structure and Energetics 

The structural characterization from an earlier high-resolution synchrotron X-ray 
diffraction^ study concluded that a-AlHs has a rhombohedral lattice of space group R3c 
(No. 167). Recent diffraction experiments^i^^ determined that 7-AIH3 has an orthorhombic 
unit cell with space group Pnnm (No. 58). The unit cell for 7-AIH3 is shown in Figure [T](a), 
compared with the hexagonal a!-AlH3 in Figure [T](c). The experimentally reported lattice 
constants are listed in Table I. The building element for both phases is the AlHg octahedron, 
where one Al atom is surrounded by six H atoms. However, the packing scheme in 7-AIH3 



3 



is more complex than that in a-AlHs. The AlHg octahedra are connected simply by sharing 
vertices in a-AlHs, as illustrated in Figure [Tt^c). The network of these octahedra produces 
only one type of Al-H-Al bridge bond which has a bond angle and a bond length of 142° and 
1.712 A, respectively. In contrast, in 7-AIH3 the AlHg octahedra are connected not only 
by sharing vertices as in a-AlHs, but also by sharing edges, as shown in Figure [T]^a). As a 
consequence, two nonequivalent Al atoms. All and A12, are created. Four nonequivalent H 
atoms (Hl-H4)^ are also identified as shown in Figure [U^a), while only one type H atom 
exists in a-AlHa. All is involved with the normal bridge bond with H, which is similar to 
that in a-AlHa but with a slight difference in the bond length and angle. However, A12 
involves a new type of double-bridge configuration A12-2H3-A12, as shown near the center of 
unit cell in Figure [T|^a). This double-bridge bond gives a smaller distance of 2.60 A between 
the two A12 atom compared with the Al-Al separation of 3.24 A in a-AlHs and 2.86 A in 
Al metal. In addition to the new double-bridge configuration, a hexagonal-ring structure is 
found consisting of two A12, two All, and four H4 atoms, as shown in Figure [T](b). The four 
Al atoms are on one plane parallel to the c-axis while the H4 atoms have a slight displace- 
ment (~0.14 A) out of this plane. These rings are connected to form linear chains along 
c-axis. The structures of double bridge in the ab-plane and the hexagonal-ring along c-axis 
are unique to 7-AIH3 and have not been found in any other hydrogen- containing aluminum 
compounds. Furthermore, the a phase possesses a higher order of symmetry than that of 
the 7 phase, with a smaller primitive trigonal unit cell. The total number of formula units 
(f.u.) in a unit cell is six and two for the 7 and a phases, respectively. The molecular volume 
(the density) for the 7 phase is found to be higher (lower) than that of the a phase by 11% 
(10%). 

Using the experimentally established space groups and unit cells as shown in Figure 
[H the total energies and structural parameters for both the a and 7 phases are studied 
from first-principles calculations. The calculated results for the fully relaxed structures are 
summarized in Table I. The theoretical lattice parameters are in good agreement with the 
observed values for both 7- and a-AlHs and are consistent with the results from previous a- 
AIII3 calculations.—!^ The calculated cohesive energies indicate that 7-AIH3 is energetically 
less stable than a-AlHa by 2.3 KJ/mol, which is in fair agreement with the measured value 
of 4.3 KJ/mol.-i^ In addition, in 7-AIII3 various bond lengths for non-equivalent Al and 
H atoms are evaluated and compared with recent experimental results of synchrotron X- 
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FIG. 1: (Color Online) Crystal structures of (a) orthorhombic 7-AIH3, (b) hexagonal-ring structure 
in 7-AIH3, and (c) hexagonal a-AlHa. The large and small spheres denote Al and H atoms, 
respectively. 

TABLE I: Calculated cohesive energies and structural parameters for 7-AIH3 and a-AlHs compared 
with values from experiments and previous calculations. 



Configuration Cohesive Energy Lattice Parameters 







eV/(f.u.) 


a(A) 


b(A) 


c(A) 





pn 


7r) 


a-AlHa: 


















Current work 




-14.052 


4.49 


4.49 


11.80 


90.0 


90.0 


120.0 


Previous work 


(Ref.iS) 




4.42 


4.42 


11.80 


90.0 


90.0 


120.0 




(Ref.i4) 




4.489 


4.4489 


11.820 


90.0 


90.0 


120.0 


Experiment 


(Ref.4) 




4.449 


4.449 


11.813 


90.0 


90.0 


120.0 


7-AIH3: 


















Current work 




-14.028 


5.43 


7.40 


5.79 


90.0 


90.0 


90.0 


Experiment 


(Ref.i6) 




5.3806 


7.3555 


5.77509 


90.00 


90.00 


90.00 




(Ref.i?) 




5.3672 


7.3360 


5.7562 


90.00 


90.00 


90.00 



ray powder diffraction^ and powder neutron diffractioni^ in Table IL Good agreements 
are found between the calculated and observed values. Without including the zero point 
energy, the calculated reaction enthalpies for equation (1) of a and 7 phases are 9.2 and 
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TABLE II: Interatomic distances (A) and bond angles (deg) obtained from the fully relaxed 
structure of 7- AIH3 . The values from the synchrotron X-ray diffractiori^^ and the powder neutron 
diffractior>i^ are also included for comparisons. 



Bond 


Calc. 


Ref.l6/Ref.l7 


Bond 


Calc. 


Ref.l6/Ref.l7 


Angle 


Calc. 


Ref.l6/Ref.l7 


A11-A12 


3.18 


3.1679/3.155 


A12-H1 


1.69 


1.668/1.657 


A11-H4-A12 


134.7 


124.0/133.88 


A12-A12 


2.62 


2.602/2.585 


A12-H2 


1.70 


1.664/1.678 


A11-H2-A12 


168.8 


171.0/169.99 


A11-H2 


1.76 


1.769/1.764 


A12-H3 


1.72 


1.70/1.755 


A12-H1-A12 


180.0 


180.0/179.99 


AI1-H4 


1.72 


1.784/1.696 


A12-H4 


1.72 


1.790/1.733 


AI2-H3-A12 


97.3 


100.7/97.53 








H3-H3 


2.30 


2.16/2.266 


H4-A11-H4 


85.9 


76.34/85.00 














H4-A12-H4 


176.3 


170.63/176.84 



10.8 KJ/mol-H2, respectively, which are also consistent with the measured values.-"^ 
B. Electronic Structure 

The electronic structure of the 7 - AIH3 is first analyzed by examining the charge distri- 
bution and charge transfer. The difference charge density Ap(r) is the difference between 
the total charge density of the solid and a superposition of the atomic charge with the same 
spatial coordinates as in the solid: 

Ap(r) = PsoHdir) - J2 Pitomir ~ Ri), (2) 

i 

where the sum is over all the atoms. Figure [2] shows such plots in the planes containing two 
different types of Al-H configurations in 7-AIH3. The solid squares and circles represent Al 
and H positions, respectively. The (002) plane shown in Figure EJ^a) contains one double- 
bridge bond A12-2H3-A12 and two other H atoms (HI and H2), while the (001) plane shown 
in Figure [2](b) contains a normal bridge H2-A11-H2. The difference density plot shows the 
positive values (solid contours) at the H positions, indicating a charge transfer from Al to 
H. As a result, aluminum is positively and hydrogen is negatively charged. The maximum 
of contour lines is in the order of 0.014 electrons/ (A) ^ for both (a) and (b). The minimum 
is about -0.002 electrons/ (A) ^ and the step size of the contours is 0.001 electrons/ (A) ^. The 
zero difference charge density line forms a closed contour around H, leaving a negative charge 
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density in the interstitial regions. 

In Figure [2]J^a), note that the positive contours of difference charge density are not exactly 
centered at the two H3 sites. The local maximum near each H3 is slightly shifted toward each 
other and the zero difference charge density line encloses both H3, showing some interactions 
between the H3 pair. In contrast, this is not seen for other H atoms, nor in the normal H-Al-H 
structure in a-AlHa. The separation between two H3 atoms in a double-bridge configuration 
is 2.3 A, which is slightly smaller than that of other neighboring H pairs. The distance of the 
two A12 atoms (2.62 A) is also small compared with that in Al metal (2.86 A). The contour 
lines in Figure [2]^a) also shows enhanced interactions between H-s and Al-d electrons. For 
comparison, the difference charge density for a-AlHs is also calculated and plotted in Figure 
E] for the (010) plane containing normal Al-H bonds. It is found that the normal bridge 





( a ) (002) plane 



(b) (001) plane 



FIG. 2: Difference charge density plots (see text) for 7-AIH3 in the planes containing Al and H 
atoms: (a) the (002) plane and (b) the (001) plane. The solid squares and circles represent Al and 
H positions, respectively. The density contour interval is 0.001 electrons/ (A"^). Charge deficiency 
is represented by dashed lines, while the density increase near the hydrogen atoms is represented 
by solid lines. 
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bond in 7-AIH3 [Figure [2]^b)] has similar characteristics as that in a-AlHs (Figure [3]). The 
binding between Al and H atoms involves a charge transfer in both a— and 7-AIH3. 

The total electronic densities of states (DOS) for both structures are calculated and given 
in Figure m The projected DOS onto H and Al for a-AlHs are also presented in Figure HI 




FIG. 3: Charge-density difference plot for a-AlHs in a plane containing Al and H atoms. The solid 
square and circle represent Al and H positions, respectively. The density contours in an interval of 
0.001 electrons/ (A^) are presented for the (010) plane. Charge deficiency is represented by dashed 
lines, while the density increase near the hydrogen atoms is represented by solid lines. 
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The angular-momentum projected DOS is evaluated by integrating over a sphere centered 
at each atom with a radius of 1.0 A for Al and 0.9 A for H, respectively. The choice of the 
H radius is based on the charge distribution around H in order to catch the charge transfer 
to H. The radius of Al is then chosen according to the Al-H bond length in order to cover 



0.5 n 




-5 5 



Energy (eV) 

FIG. 4: (Color Online) Calculated electronic density of states (DOS) projected onto (a) H and (b) 
Al in a-AlHs, as well as the total DOS for the (c) a-AlHs and (d) 7-AIH3. The Fermi level is at 
zero. 
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as much of the interstitial region as possible without overlapping the spheres. The Fermi 
energy is set to zero in these plots. 

The current results for a-AlHs are comparable with previous calculations.— >i^>i^ The H-s 
component spans the whole energy range. The lowest peak from ~ -8.7 to -4.0 eV is found 
to correspond to H-s and Al-s states, while the higher-energy features are composed of H-s 
and Al-p or -d states. The total bond width is similar in both a- and 7-AIH3. A difference 
between the two is clearly seen near -7 eV, where a small band gap is created in the low 




FIG. 5: (Color Online) Calculated electronic density of states projected onto non-equivalent atoms 
in 7-AIH3. The Fermi energy is set to be zero. 
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energy region in 7-AIH3. This new gap gives rise to a broad, separated DOS peak from -8.7 
to -7.0 eV. 

In order to understand the new feature found in the DOS of 7-AIH3, the projected density 
of states onto non-equivalent Al and H atoms are analyzed in Figure O The same radius 
values used in a-AlHs are adopted for integrating over the spheres around Al and H atoms. 




(002) plane 

FIG. 6: Total charge density for the lowest two energy bands projected onto the (002) plane, 
the same plane as shown in Figure [2]J^a) . The solid circles represent H positions. The maximum, 
minimum and the interval of the contours are 0.85, 0.002, and 0.1 electrons/ (A^), respectively. 
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The valence states below the Fermi energy can be grouped into three energy regions: the 
low( -8.7 to -7.0), the middle(-6.8 to -4.0), and the high( -4.0 to eV) energy regions. 
The broad peak in the low energy region of -8.7 to -7.0 is composed of considerable 




FIG. 7: (Color Online) Total charge density isosurface for the lowest two energy bands in 7-AIH3. 
The cell structure is the same as shown in Figure [T]^a). The constant charge density surface shown 
is for 0.45 electrons/(A)^. The charge accumulation within the hexagonal-ring structure defined in 
Figure [T][b) is clearly shown. 
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H3-S states, with some s contributions from other H and two Al atoms. This broad peak 
consists of the first two energy bands separated from the other bands at higher energies by 
a small gap. In Figure El the total charge density for these two bands is calculated and 
plotted for the (002) plane, the same plane shown in Figure [2|^a). The maximum, minimum 
and interval of the contours are 0.85, 0.002, and 0.1 electrons/ (A) ^, respectively. It shows 
noticeable charge accumulation around the H3-H3 pair. The three-dimensional isosurface 
of the total charge density for the same lowest two bands is shown in Figure U\ which has 
a density value of 0.45 electrons/ (A) ^. Again, the interaction between the H3-H3 pair is 
illustrated. In addition, considerable charge is found around four H4 atoms within the same 
hexagonal-ring as defined in Figure [T](b). Therefore, the broad peak in the low energy region 
is attributed mainly to the unique double-bridge and hexagonal-ring structures in 7-AIH3. 

In fact, the projected DOS associated with H3 is reduced in the middle energy region 
(-6.8 to -4.0 eV) compared with other H atoms. The missing amplitude is shifted to both 
the low (-8.7 to -7.0 eV) and high (-4.0 to eV) energy regions. The former contains 
an increased H3-H3 interaction, the H-s and A12-s coupling, while the latter exhibits an 
enhanced interaction between H3-s and A12-p and -d states. The stronger interaction of 
H3-S states with higher-energy A12-p and -d states due to geometry gives rises to a higher 
electronic energy for the system. 

C. Vibrational Properties 

It is expected that the high vibrational frequencies of light hydrogen atoms may result 
in a significant zero-point energy, which needs to be included in studying the ground-state 
properties. In order to consider the zero-point energy E^pt contributions to the energetics 
between the two AIH3 phases, the phonon density of states (DOS) are evaluated for both 
a- and 7-AIII3. Although the absolute value of E^pt is noticeable for the two phases, the 
difference turns out to be small AE^pt = Ezpt{oi) — Ezpti'j) = 0.2 KJ/mol. By including this 
correction to the total energies previously listed in Table I, the 7 phase has a slightly higher 
energy than the a phase by 2.1 KJ/mol. 

The partial phonon density of states (DOS) for atom r is defined as: 

37V 
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FIG. 8: (Color Online) Total and partial phonon density of states (DOS) for a-AlHs. 
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where the N is the total number of atoms per unit cell, q is the phonon momentum, j labels 
the phonon branch, e^(q, j) is the phonon polarization vector for atom r, and a;(q, j) is the 
phonon frequency. The total and projected phonon density of states for a- and 7- AIH3 are 
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FIG. 9: (Color Online) Total and partial phonon density of states (DOS) for 7-AIH3 
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presented in the Figures [8] and [H respectively. Our calculated phonon DOS for a-AlHs is in 
good agreement with previous calculations,— showing three distinct groups in the whole 
frequency range. The decomposed phonon DOS indicates that the low frequency modes 
below 350 cm"^ are mainly from Al atoms, while the middle and high frequency modes 
are H- dominated. This is due to the large mass ratio between H atom to Al atom. The 
high-frequency phonons above 1550 cm~^ are associated with H motion in the Al-H bond 
stretching modes. No vibrational modes was found in the frequency region from 1025 to 
1550 cm~^, yielding gap of 525 cm~^ in a-AlHa. 

The phonon DOS for 7-AIH3 is more complex than that for a-AlHa due to the structural 
difference. The 7- phase unit cell has much more atoms (24) than the a-phase (8), therefore 
the calculated phonon DOS for 7-AIH3 in Figure exhibit more features. Comparing with 
the phonon DOS of a-AlHs, although the modes located in the three frequency regions 
remain, multiple new additional peaks appear in the spectrum. In order to understand how 
different atomic species contribute to these vibrations, the partial phonon DOS is calculated 
for the four types of H and two types of Al atoms, as illustrated in Figure [9l 

The new features labeled in Figure [9] can be summarized as follows: (1) four new vibra- 
tional peaks appear within the gap region of the C1-AIH3 spectrum in the frequency range 
of 1200 to 1500 cm~^. These are vibrational modes of H3 in the plane of the A12-2H3-A12 
complex. It consists of both bond stretching and bending motions by 113 on this plane. (2) 
A new strong narrow peak is introduced in the middle frequency region around a frequency 
of 875 cm~^, corresponding to the H3 vibrations perpendicular to the A12-2H3-A12 plane. 
(3) An additional broaden peak appears in the low frequency region from 375 to 475 cm~^, 
which is dominated by A12 in-plane vibrational modes, coupled with some contributions 
from HI and H2 atoms with little contribution from H3. The displacements of the two 
A12 atoms have the same magnitude but opposite directions, indicating a paired motion 
between them. As a consequence, these modes have higher frequencies than the other Al 
modes. The first three new features are largely related to the motions of the two H3 and the 
two A12 atoms in the double-bridge structure. (4) A set of two peaks around 1550 cm~^ are 
associated with the out-of-plane motion of four 114 atoms in the hexagonal-ring structure 
shown in Figure [T](b). (5) The set of high-frequency peaks around 2012 cm~^ are associated 
with HI atoms which connect two A12 atoms in opposite directions. The angle A12-H1-A12 
is 180°. The vibrational modes of HI are in the ab plane and along the A12-H1-A12 line. 
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These new features displayed in the 7-AIH3 phonon DOS can serve as indicators of these 
special Al and H arrangements. 



IV. CONCLUSION 



We have performed pseudopotential density-functional calculations to study the elec- 
tronic and phonon vibrational properties for the newly reported aluminum hydride 7-AIH3 
structure. The calculated structural parameters are in good agreement with results from 
the diffraction experiments. Our energetic study for the AIH3 system indicates that 7-AIH3 
is less stable than a;-AlH3 by ~ 2 KJ/mol. The unique double-bridge configuration in 7- 
AIH3 is investigated by examining the electronic properties and phonon vibrational modes. 
It was found that the double-bridge arrangement modified the binding between H and Al 
atoms. The projected DOS indicates that interactions exist between the double-bridge H3 
atoms and that the interaction between H3-s and higher-energy A12-p and -d states is en- 
hanced. The latter yields a higher electronic energy in 7-AIH3. New features in phonon 
vibrational spectrum associated with double-bridge bonds and hexagonal-ring complex were 
also identified. 
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